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Abstract 

We apply the thermodynamic Bethe Ansatz to investigate the high energy 
behaviour of a class of scattering matrices which have recently been proposed 
to describe the Homogeneous sine-Gordon models related to simply laced Lie 
algebras. A characteristic feature is that some elements of the suggested S- 
matrices are not parity invariant and contain resonance shifts which allow 
for the formation of unstable bound states. From the Lagrangian point of 
view these models may be viewed as integrable perturbations of WZNW-coset 
models and in our analysis we recover indeed in the deep ultraviolet regime 
the effective central charge related to these cosets, supporting therefore the 
S-matrix proposal. For the S'?7(3)/ c -model we present a detailed numerical 
analysis of the scaling function which exhibits the well known staircase pat- 
tern for theories involving resonance parameters, indicating the energy scales 
of stable and unstable particles. We demonstrate that, as a consequence of 
the interplay between the mass scale and the resonance parameter, the ul- 
traviolet limit of the HSG-model may be viewed alternatively as a massless 
ultraviolet-infrared- flow between different conformal cosets. For k = 2 we 
recover as a subsystem the flow between the tricritical Ising and the Ising 
model. 
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1 Introduction 

The thermodynamic Bethe ansatz (TBA) is established as an important method 
which serves to investigate "off-shell" properties of 1+1 dimensional quantum field 
theories. Originally formulated in the context of the non-relativistic Bose gas by 
Yang and Yang Q, it was extended thereafter by Zamolodchikov |2| to relativistic 
quantum field theories whose scattering matrices factorize into two-particle ones. 
The latter property is always guaranteed when the quantum field theory in question 
is integrable. Provided the S-matrix has been determined in some way, for instance 
via the bootstrap program || or by extrapolating semi-classical results, the TBA 
allows to calculate the ground state energy of the integrable model on an infinite 
cylinder whose circumference is identified as compactified space direction. When 
the circumference is sent to zero the effective central charge of the conformal field 
theory (CFT) governing the short distance behaviour can be extracted. In the case 
in which the massive integrable field theory is obtained from a conformal model 
by adding a perturbative term which breaks the conformal symmetry, the TBA 
constitutes therefore an important consistency check for the S-matrix. 

The main purpose of this manuscript is to apply this technique to a class of scat- 
tering matrices which have recently been proposed Q to describe the Homogeneous 
sine-Gordon models (HSG) [(J related to simply laced Lie algebras. The latter 
have been constructed as integrable perturbations of WZNW-coset theories of 
the form G^/H, where G is a compact simple Lie group, H C G a maximal abelian 
torus and k > 1 an integer called the "level". These models constitute particular 
deformations of coset-models @, where the specific choice of the groups ensures 
that these theories possess a mass gap ||. The defining action of the HSG-models 
reads 



Here Soft denotes the coset action, ( , ) the Killing form of G and g(x) a group 
valued bosonic scalar field. A± are arbitrary semi-simple elements of the Cartan 
subalgebra associated with H, which have to be chosen not orthogonal to any 
root of G and play the role of continuous vector coupling constants. The latter 
constraints do not restrict the parameter choice in the quantum case with regard to 
the proposed S-matrix which makes sense for every choice of A±. They determine 
the mass ratios of the particle spectrum as well as the behaviour of the model under 
a parity transformation. The parameters m and f3 2 = 1/k + 0(l/k 2 ) are the bare 
mass scale and the coupling constant, respectively. The non-perturbative definition 
of the theory is achieved by identifying (A + , g(x)~ 1 A_g(x)) with a matrix element 
of the WZNW-field g(x) taken in the adjoint representation, which is a spinless 
primary field of the coset-CFT and in addition exchanging (3 by 1/k and m by the 
renormalised mass ||. Some of the conformal data of Scft[<?]) which are in principle 
extractable from the TBA analysis are the Virasoro central charge c of the coset 
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and the conformal dimensions A, A of the perturbing operator in the massless limit 



k dim G „ k — 1 , „ „ t h 

* k+h k + h k+h w 

Here I denotes the rank of G and h its Coxeter number. Since we have A < 
1 for all allowed values of k, the perturbation is always relevant in the sense of 
renormalisationQ. 

The simplest example of a HSG theory is the complex sine-Gordon model || |J 
associated with the coset SU{2)k/U{l). As we will argue below, more complicated 
HSG theories can be viewed as interacting copies of complex sine-Gordon theo- 
ries. The classical equations of motion of these models correspond to non-abelian 
affine Toda equations || [Kj], which are known to be classically integrable and ad- 
mit soliton solutions. Identifying these solutions by a Noether charge allows for a 
semi-classical approach to the quantum theory by applying the Bohr-Sommerfeld 
quantization rule. The integrability on the quantum level was established in || by 
the construction of non-trivial conserved charges, which suggests the factorization 
of the scattering matrix. Based on the assumption that the semi-classical spectrum 
is exact, the S-matrix elements have then been determined in Q by means of the 
bootstrap program for HSG-models related to simply laced Lie algebras. 

The proposed scattering matrix consists partially of i copies of minimal su(k)- 
affine Toda field theories (ATFT) whose mass scales are free parameters. The 
scattering between solitons belonging to different copies is described by an S-matrix 
which violates parity ||. These matrices possess resonance poles and the related 
resonance parameters which characterize the formation of unstable bound states 
are up to free choice. In the TBA-analysis these resonances lead to the "staircase 
patterns" in the scaling function, which have been observed previously for similar 
models 1T3]. However, in comparison with the models studied so far, the HSG 



models are distinguished in two aspects. First they break parity invariance and 
second some of the resonance poles can be associated directly to unstable particles 
via a classical Lagrangian. 

One of the main outcomes of our TBA-analysis is that the suggested [[|] scatter- 
ing matrix leads indeed to the coset central charge (0), which gives strong support 
to the proposal. 

In addition, we present a detailed numerical analysis for the >S?7(3)-HSG model, 
but expect that many of our findings for that case are generalizable to other Lie 
groups. The presence of two parameters, i.e. the mass scale and the resonance 
parameter allow, similar as for staircase models studied previously, to describe 
the ultraviolet limit of the HSG-model alternatively as the flow between different 
conformal field theories in the ultraviolet and infrared regime. We find the following 



*We slightly abuse here the notation and use CG k instead of c Gk /u^iy. Since we always en- 
counter these type of coset in our discussion, we can avoid bulky expressions in this way. 
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massless flow 



UV = SU(3) k /U(l) 2 <-> SU{2) k /U{l) ® SU{2) k /U{l) = IR . (3) 

We also observe the flow (SU(3) k /U(l) 2 )/(SU(2) k /U(l)) -> SU{2) k /U{\) as a 
subsystem inside the HSG-model. For k = 2 this subsystem describes the flow 



between the tricritical Ising and the Ising model previously studied in |Tj|. In 
terms of the HSG-model we obtain the following physical picture: The resonance 
parameter characterizes the mass scale of the unstable particles. Approaching the 
extreme ultraviolet regime from the infrared we pass various regions: At first all 
solitons are too heavy to contribute to the off-critical central charge, then the two 
copies of the minimal ATFT will set in, leading to the central charge corresponding 
to IR in (|3|) and finally the unstable bound states will start to contribute such that 
we indeed obtain (Q) as the ultraviolet central charge of the HSG-model. 

The two values of the resonance parameter and oo are special, corresponding 
in the classical theory to a choice of the vector couplings in (|I]) parallel to each other 
or orthogonal to a simple root, respectively. In the former case parity is restored on 
the classical as well as on the TBA-level and the central charge corresponding to 
UV in (^j is also recovered, whereas in the latter case the two copies of the minimal 
ATFT are decoupled and unstable bound states may not be produced leading to 
the central charge IR in ([3]). 

Our manuscript is organized as follows: In section 2 we briefly recall the main 
features of the two-particle HSG S-matrix elements stating them also newly in 
form of an integral representation. In particular, we comment on the link between 
unstable particles and resonance poles as well as on the loss of parity invariance. In 
section 3 we introduce the TBA equations for a parity violating system and carry 
out the ultraviolet limit recovering the expected coset central charge. In section 4 
we present a detailed study for the SU(3) k — HSG model. We discuss the staircase 
pattern of the scaling function and illustrate how the UV limit for the HSG-model 
may be viewed as the UV-IR flow between different conformal models. We extract 
the ultraviolet central charges of the HSG-models. We study separately the case 
when parity is restored, derive universal TBA-equations and Y-systems. In section 
5 we present explicit examples for the specific values k = 2,3, 4, oo. Our conclusions 
are stated in section 6. 



2 The homogeneous sine-Gordon S-matrix 

We shall now briefly recall the main features of the proposed HSG scattering matrix 
in a form most suitable for our discussion. Labelling the solitons by two quantum 
numbers, we take the two-particle scattering matrix between soliton (a, i) and soli- 
ton (b,j), with 1 < a,b < k — 1 and 1 < i, j < £, as a function of the rapidity 
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difference 9 to be of the general form S l J b (8). The particular structure of the con- 
jectured HSG S-matrix makes it suggestive to refer to the lower indices as main 
quantum numbers and to the upper ones as colour. In ||] it was proposed to de- 
scribe the scattering of solitons which possess the same colour by the S-matrix of 
the Zfc-Ising model or equivalently the minimal su(k)-ATFT [ |TT|j 

min(a,fe) — 1 

S5,{0) = (a + b) e (\a-b\)e J] (a + b-2n) 2 e (4) 

71=1 

-1 



/dt 7lt ( lit \ -,n , , 

— 2 cosh— ^2 cosh— -I J e~ lW . (5) 



ab 

Here we have introduced the abbreviation (x)g = sinh ^{d + i^rr)/ s ^ nn |(^ — *lf) ^ or 
the general building blocks and denote the incidence matrix of the sii(/c)-Dynkin 
diagram by I. We re-wrote the above S-matrix from the block form into a form 
of an integral representation ([5]), since the latter is more convenient with respect to 
the TBA analysis. This calculation may be performed by specializing an analysis 
in O, [15| to the particular case at hand. The scattering of solitons with different 



colour quantum numbers was proposed to be described by 

min(a,6)— 1 

Slid) = (V ij ) ah II Ha-b\-l-2n) 6+rTij , 0,2 (6) 



n=0 

ab 



(%-rexp 



J j (2 cosh j - /) 1 e- u ^\ Kl ± 0, 2 , (7) 



with i^ 9 denoting the Cartan matrix of the simply laced Lie algebra g. Here the 
Vij = rfi are ar bitrary k-th roots of —1 taken to the power a times b and the shifts 
in the rapidity variables are functions of the vector couplings cr^, which are anti- 
symmetric in the colour values = —<Jji- Due to the fact that these shifts are 
real, the function S^ b (6) for i 7^ j will have poles beyond the imaginary axis such 
that the parameters characterize resonance poles. An important feature is that 
is not parity invariant, where parity is broken by the phase factors r] as well as 
the shifts o . As a consequence, the usual relations 

sm = sue) = s: b (-e*y and s£(*)5£(-0) = 1 (8) 

satisfied by the parity invariant objects are replaced by 

s l M = sa-er and s%{e)st{-e) = 1 (9) 

for the scattering between solitons with different colour values. Important to note 
is that the first equality in (§) has no analogue in @. Thus, instead of being real 
analytic, as S™ b (0), the parity violation forces Hermitian analyticity of S l J b (9) for 
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% ^ j. Anti-particles are constructed in analogy to the ATFT, that is from the 
automorphism which leaves the sw(fc)-Dynkin diagram invariant, such that (a, i) = 
(k — a,i). The colour of a particle and its anti-particle is identical. The crossing 
relation of the S-matrix then reads 

Sm = 4- a)b (0) = SHiin - 6) . (10) 



For a general and more detailed discussion of these analyticity issues see |16| and 
references therein. 

Analyzing the above S-matrix we have the following picture concerning the 
formation of bound states: Two solitons with the same colour value may form a 
bound state of the same colour, whilst solitons of different colour with Ky ^ 0, 2, 
say (a, i) and (6, j), may only form an unstable state, say (c, k) whose lifetime and 
energy scale are characterized by the parameter o by means of the Breit-Wigner 



formula, see e.g. [17], in the form 



(T k ) 2 

(Mf) 2 -^- = (M t a ) 2 + (Ml) 2 + 2M l a M J b coshacose (11) 
M|r| = 2M*M b i sinh|a|sin6 , (12) 

where the resonance pole in S^ b (9) is situated at Or = a — iO and T~ denotes 
the decay width of the unstable particle with mass Mf . In the case a = b these 



unstable states can be identified with solitons in the semi-classical limit ||, |24 
When a becomes zero, (|r2|) shows that the unstable particles become stable, but 
are still not at the same footing as the other asymptotically stable particles. They 
become virtual states characterized by poles on the imaginary axis beyond the 
physical sheet. 

How many free parameters do we have in our model? Computing mass shifts 
from renormalisation, we only accumulate contributions from intermediate states 
having the same colour as the two scattering solitons. Thus, making use of the well 
known fact that the masses of the minimal sw(/c)-affine Toda theory all renormalise 



with an overall factor Ul8| , i.e. for the solitons (a, i) we have that 5M l a /M l a equals a 



constant for fixed colour value i and all possible values of the main quantum number 
a, we acquire in principle I different mass scales mi, . . . ,m£ in the HSG- model. In 
addition there are £ — 1 independent parameters in the model in form of the possible 
phase shifts = for each i,j such that Kf- ^ 0,2. This means overall we 
have 2^—1 independent parameters in the quantum theory. There is a precise 
correspondence to the free parameters which one obtains from the classical point of 
view. In the latter case we have the 2£ independent components of A± at our free 
disposal. This number is reduced by 1 as a result of the symmetry A + — > cA + and 
A_ — > c _1 A_ which introduces an additional dependence as may be seen from the 
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explicit expressions for the classical mass ratios and the classical resonance shifts 



Mi _ \{a i -K + ){a i - A_) _ ^ / (a, ■ A + )(^ ■ A_) ^ 



Mi V ( a i " A ~) ( a i " A +) ' K V ("i • A -) ( a S ■ A +) 



Here the on are simple roots. 

In comparison with other factorizable scattering matrices involving resonance 
shifts, studied in the literature so far, the proposed HSG scattering matrices differ 
in two aspects. First of all, they are not parity invariant and second they allow to 
associate a concrete Lagrangian description. The latter fact can be used to support 
the picture outlined for the full quantum field theory by a semi-classical analysis. 



In 24 the semi-classical mass for the soliton fa, i) was found to be 



rrij na 

M = — y sin — 14 

tt/3 2 k 



where j3 is a coupling constant and the rrii are the £ different mass scales. 



3 TBA with parity violation and resonances 

In this section we are going to determine the conformal field theory which governs 
the UV regime of the quantum field theory associated with the S-matrix elements 
(fD and (H) . According to the defining relation ([I]) and the discussion of the previous 
section, we expect to recover the WZNW-coset theory with effective central charge 
(0) in the extreme ultraviolet limit. It is a well established fact that such high 
energy limits can be performed by means of the TBA. Since up to now such an 
analysis has only been carried out for parity invariant S-matrices, a few comments 
are due to implement parity violation. The starting point in the derivation of the 
key equations are the Bethe ansatz equations, which are the outcome of dragging 
one soliton, say of type A = (a,i), along the world line. For the time being we do 
not need the distinction between the two quantum numbers. On this trip the formal 
wave-function of A picks up the corresponding S-matrix element as a phase factor 
when meeting another soliton. Due to the parity violation it matters, whether the 
soliton is moved clockwise or counter-clockwise along the world line, such that we 
end up with two different sets of Bethe Ansatz equations 

e iLM AS mte A "Q Sab (q a _ 0b) = 1 and e -iLM A sin^ A "Q g^^ - Q A ) = I , 

(15) 

with L denoting the length of the compactified space direction. These two sets of 
equations are of course not entirely independent and may be obtained from each 
other by complex conjugation with the help of relation (||]). We may carry out the 
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thermodynamic limit of ( |I5D in the usual fashion |2| and obtain the following sets 
of non-linear integral equations 

t\( e ) + J2 ®ab*L+(6) = r Ma cosh 9 (16) 

B 

£a(8)+J2® ba * L b( 9 ) = rM A cosh6 . (17) 

B 

As usual we let the symbol '*' denote the rapidity convolution of two functions 
defined by / * g{9) := J d9' /2ix f(9 — 9')g{9'). Here r = m\T~ x is the inverse 
temperature times the overall mass scale m\ of the lightest particle and we also re- 
defined the masses by M l a —>■ M l a jm\ keeping, however, the same notation. As very 
common in these considerations we also introduced the so-called pseudo-energies 
6 a(@) = e A.(~@) an d the related functions L A (9) = ln(l + e -e ^w). The kernels in 
the integrals carry the information of the dynamical interaction of the system and 
are given by 

$ab(0) = $ba(-0) = -i-^lnS A B{e) ■ (18) 

dv 

The statistical interaction is chosen to be of fermionic type. Notice that (O) may 
be obtained from ([TJ]) simply by the parity transformation 9 — > —9 and the first 
equality in fll8D . The main difference of these equations in comparison with the 
parity invariant case is that we have lost the usual symmetry of the pseudo-energies 
as a function of the rapidities, such that we have now in general e\{9) ^ e^(6 l ). 
This symmetry may be recovered by restoring parity. 

The scaling function, which can also be interpreted as off-critical Casimir energy, 
may be computed similar as in the usual way 

CO 

c(r) = ^^M A J d9 coshe (LJ(9) + L+(6)) , (19) 
A o 

once the equations fllSD have been solved for the e A (9). Of special interest is the 
deep UV limit, i.e. r — > 0, of this function since then its value coincides with the 
effective central charge c e fr = c — 12(A + A ) of the conformal model governing the 
high energy behaviour. Here c is the Virasoro central charge and A , A are the 
lowest conformal dimensions related to the two chiral sectors of the model. Since the 
WZNW-coset is unitary, we expect that A , A = and c c s = c. This assumption 
will turn out to be consistent with the analytical and numerical results. 

3.1 Ultraviolet central charge for the HSG model 

In this section we follow the usual argumentation of the TBA-analysis which leads 
to the effective central charge, paying, however, attention to the parity violation. 
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We will recover indeed the value in (0) as the central charge of the HSG-models. 
First of all we take the limits r, 9 — > of ( |16|) and (|T7|) . When we assume that 



the kernels $ab(#) are strongly peaked[] at 9 = and develop the characteristic 
plateaus one observes for the scaling models, we can take out the L-functions from 
the integral in the equations ([161), (|T7|) and obtain similar to the usual situation 

oo 

e±(0) + V N AB L±(0) = with N AB = ±- [ d9 <b AB (9) . (20) 

b 271 J 

— oo 

Having the resonance parameter a present in our theory we may also encounter 
a situation in which & AB (9) is peaked at 9 = ±o~. This means in order for (|20|) 
to be valid, we have to assume e^(0) = e^(±cr) in the limit r — > in addition 
to accommodate that situation. For the last assumption we will not provide a 
rigorous analytical argument, but will justify it instead for particular cases from 
the numerical results (see e.g. figure 1). Note that in ( p0|) we have recovered the 
parity invar iance. 

For small values of r we may approximate, in analogy to the parity invariant sit- 
uation, rM^cosh^ by r /2M^exp 9, such that taking the derivative of the relations 
(|T6|) and fllT| ) thereafter yields 



oo 

m + _L V f d9 > ^l e l d ff^ * -Ma exp 9 . (21) 



B -oo 



l + exp(eg(0')) ^ 
The scaling function acquires the form 

oo 

c(r) ~ |^ V M4 /" d0 exp (L^(0) + L+(0)) , for r ~ (22) 
2tt z z — ' / 



in this approximation. Replacing in ( p2|) the term r/2M^exp0 by the l.h.s. of (|21| ) 
a few manipulations lead to 



hmc(r)~ — J2J2 de A 

^ — I A J 



£(«>) 



A 4(0) 



ln(l + exp(-e^)) + 



A 



1 + exp(e 



(23) 



TThat this assuption holds for the case at hand is most easily seen by noting that the logarithmic 
derivative of a basic building block (x)$ of the S-matrix reads 



d9 



Sill T X ^-^ / 7T 

-i— ]n(x) e = r = -2 > s in I - x e 

V ' cosh6»-cosf^x) ^ \k 



-n\8\ 



n=l 

From this we can read off directly the decay properties. 
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Upon the substitution y p A = 1/(1 + exp(e^)) we obtain the well known expression 
for the effective central charge 



^E< 1+e 4s (0)) ) ■ (*) 

Here we used the integral representation for Roger's dilogarithm function C(x) = 

V 2 Jo dvi^vliv - l ) - M 1 - v)/v), and the facts that = e A(°)' Va(.°°) = 

y^(oo) = 0. This means we end up precisely with the same situation as in the parity 
invariant case: Determining at first the phases of the scattering matrices we have 
to solve the constant TBA-equation ( P0| ) and may compute the effective central 
charge in terms of Roger's dilogarithm thereafter. Notice that in the ultraviolet 
limit we have recovered the parity invariance and (24) holds for all finite values of 
the resonance parameter. 

For the case at hand we read off from the integral representation of the scattering 
matrices 

N^ b = 5 t3 5 ab -mK^)-\ b . (25) 

With N% in the form (g|) and the identifications Q\ = Ub^li 1 + exp(-e l b (0)))^ 1 
the constant TBA-equations are precisely the equations which occurred before in 
the context of the restricted solid-on-solid models [T^, |2(|. It was noted in there 



that (|20D may be solved elegantly in terms of Weyl-characters and the reported 
effective central charge coincides indeed with the one for the HSG-models 

It should be noted that we understand the iV-matrix here as defined in (|20| ) 
and not as the difference between the phases of the S-matrix. In the latter case we 
encounter contributions from the non-trivial constant phase factors 77. Also in that 
case we may arrive at the same answer by compensating them with a choice of a 
non-standard statistical interaction as outlined in Wl . 



We would like to close this section with a comment which links our analysis to 
structures which may be observed directly inside the conformal field theory. When 



one carries out a saddle point analysis, see e.g. [|2T|], on a Virasoro character of the 
general form 

^ (q)i ■ ■ ■ (q)( k -i)i 

with (q) m = rifcLi(l ~ q k )i one recovers the set of coupled equations as (|20| ) and the 



corresponding effective central charge is expressible as a sum of Roger's dilogarithms 
as fl2~4l) . Note that when we choose g = A\ the HSG-model reduces to the minimal 
ATFT and ( p6|) reduces to the character formulae in ]23| . Thus the equations (|20| ) 



and (24) constitute an interface between massive and massless theories, since they 
may be obtained on one hand in the ultraviolet limit from a massive model and on 
the other hand from a limit inside the conformal field theory. This means we can 
guess a new form of the coset character, by substituting (|25| ) into (f2~6l). However, 
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since the specific form of the vector B does not enter in this analysis (it distinguishes 
the different highest weight representations) more work needs to be done in order 
to make this more than a mere conjecture. We leave this for future investigations. 



4 The SU(3)a;-HSG model 

We shall now focus our discussion on G = SU(3)k- First of all we need to establish 
how many free parameters we have at our disposal in this case. According to the 
discussion in section 2 we can tune the resonance parameter and the mass ratio 

a := (72i = — o"i2 and mi/m 2 . (27) 

It will also be useful to exploit a symmetry present in the TBA-equations re- 
lated to SU(3)k by noting that the parity transformed equations (|17D turn into the 
equations (ITH) when we exchange the masses of the different types of solitons. This 
means the system remains invariant under the simultaneous transformations 

9 — > —9 and mi/m2 m 2 /m 1 . (28) 

For the special case m 1 /m 2 = 1 we deduce therefore that ei 1 ^) = ei 2 ^— 9), meaning 
that a parity transformation then amounts to an interchange of the colours. Fur- 
thermore, we see from (|17| ) and the defining relation o = oi\ = —<Ji2 that changing 
the sign of the rapidity variable is equivalent to a —>■ —a. Therefore, we can restrict 
ourselves to the choice a > without loss of generality. 



4.1 Staircase behaviour of the scaling function 



We will now come to the evaluation of the scaling function (19) for finite and 
small scale parameter r. To do this we have to solve first the TBA equations (|16|) 
for the pseudo-energies, which up to now has not been achieved analytically for 
systems with a non-trivial dynamical interaction due to the non-linear nature of 
the integral equations. Nonetheless, numerically this problem can be controlled 
relatively well. In the UV regime for r 1 one is in a better position and can 
set up approximate TBA equations involving formally massless particles^ for which 
various approximation schemes have been developed which depend on the general 
form of the L-functions. Since the latter is not known a priori, one may justify 
ones assumptions in retrospect by referring to the numerics. In section 5 we present 



■''The concept of massless scattering has been introduced originally in [jT3| as follows: The on- 
shell energy of a right and left moving particle is given by E± = M/2e ±e which is formally obtained 
from the on-shell energy of a massive particle E = mcoshfl by the replacement 9 — > 9 ± a/2 and 
taking the limit m — > 0, a —* oo while keeping the expression M = me e+a ^ 2 finite. It is these 
on-shell energies we will encounter in our analysis. 
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numerical solutions for the equations (IB ) for various levels k showing that the L 



functions develop at most two (three if m,\ and m2 are very different) plateaus in the 
range hi ~ < 9 < In - and then fall off rapidly (see figure 1). This type of behaviour 
is similar to the one known from minimal ATFT 0, |29j , and we can therefore adopt 



various arguments presented in that context. The main difficulty we have to deal 
with here is to find the appropriate "massless" TBA equations accommodating the 
dependence of the TBA equations on the resonance shifts a. 
We start by separating the kernel (|T%D into two parts 



<P ab {0) = $£(0) = Jdt [d ab - 2 cosh § (2 cosh f - e~ ite , (29) 

i>aM = *&(0 + *i<) = Jdt (2 cosh f -J);; e- ite (30) 

Here <f) a b{9) is just the TBA kernel of the su(k) -minimal ATFT and in the remaining 
kernels ip a b(9) we have removed the resonance shift. Note that <fr, ip do not depend on 
the colour values i, j and may therefore be dropped all together in the notation. The 
integral representations for these kernels are obtained easily from the expressions 
(|5|) and ([/]). They are generically valid for all values of the level k. The convolution 
term in (|T6|) in terms of 0, ip is then re- written as 



i fc-i fc-i i fc-i 

EE^* = E * L >w + E E * l w - ^) • ( 31 ) 

j=l 6=1 6=1 j=l 6=1 

These equations illustrate that whenever we are in a regime in which the second term 
in ( pTD is negligible we are left with I non-interacting copies of the su( k ) -minimal 
ATFT. 

We will now specialize the discussion on the sw(3)fc-case for which we can elim- 
inate the dependence on a in the second convolution term by performing the shifts 
9 — » 9 ± a/2 in the TBA equations. In the UV limit r —>■ with a ^> 1 the shifted 
functions can be approximated by the solutions of the following sets of integral 
equations 

fc-i fc-i 



^w+E^*^w + E^* L ?w = r ' M « e±e ( 32 ) 

6=1 6=1 

fc-1 

£(*)+X>.6 = r'il^e ±e , (33) 



6=1 



where we have introduced the parameter r' = res/2 familiar from the discussion of 
massless scattering and the masses Ma = Ma . The relationship between the 
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solutions of the massless system (|32"D, ( |33"D and those of the original TBA-equations 
is given by 

e CO/(3)(0) = £ +/-(0 T(7 /2) for ln(r/2) < ±9 < ln(r/2) + a (34) 
e ( 1 )/( 2 )(0) = r /+ (^±a/2) for ±0 « min[ln(2/r),ln(r/2) +<r)] (35) 

where we have assumed mi = m-i. Similar equations may be written down for 
the generic case. To derive ( |35|) we have neglected here the convolution terms 
(ipab * Lfr)(9 + cr) and {ip ab * L^)(6 — a) which appear in the TBA-equations for 
£a\0) and 6a (9), respectively. This is justified by the following argument: For 
\9\ 3> ln2/r the free on-shell energy term is dominant in the TBA equations, i.e. 
e l a (9) « rM" cosh 9 and the functions L l a (9) are almost zero. The kernels ip ab are 
centered in a region around the origin # = outside of which they exponentially 
decrease, see footnote in section 3.1. for this. This means that the convolution 
terms in question can be neglected safely if 9 <C ln(r/2) + a and 9 ^> ln(2/r) — a, 
respectively. Note that the massless system provides a solution for the whole range 
of 9 for non-vanishing L-function only if the ranges of validity in ( p4|) and ( |35| ) 
overlap, i.e. if ln(r/2) <C min[ln(2/r), ln(r/2) + a] which is always true in the 
limit r — > when a 0. The solution is uniquely defined in the overlapping region. 
These observations are confirmed by our numerical analysis below. 

The resulting equations ( |3"3"| ) are therefore decoupled and we can determine L + 
and L~ individually. In contrast, the equations (|3~2"D for are still coupled to each 
other due to the presence of the resonance shift. Formally, the on-shell energies for 
massive particles have been replaced by on-shell energies for massless particles in 
the sense of |13| , such that if we interpret r' as an independent new scale parameter 
the sets of equations fl32| ) and ( |33"D could be identified as massless TBA systems in 
their own right. 

Introducing then the scaling function associated with the system ( |3"2] ) as 

c °( r ') = %Y,l de K ^ L+ a( ) + M a e- e L-(9)} (36) 
71 a=i ^ 

and analogously the scaling function associated with ( [33D as 

fc— l 

J2jd9 [M+ e 9 L+(9) + M~ e' s L-(9)] (37) 



3, k 1 
r 



a=l 



we can express the scaling function fllTf ) of the HSG model in the regime r 
0, a ^> 1 approximately by 

c(r,a) = ^EE M » [ d9 [e K(9-a/2) + e- e Li(9 + a/2 

i=l,2 a=l J 

~ c (/) + c Q (/) . (3J 
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Thus, we have formally decomposed the massive S77(3)fc-HSG model in the UV 
regime into two massless TBA systems (|32|) and (|33|) , reducing therefore the problem 
of calculating the scaling function of the HSG model in the UV limit r —>■ to the 
problem of evaluating the scaling functions ( |36D and ([37]) for the scale parameter 
r'. The latter depends on the relative size of ln(2/r) and the resonance shift a/ 2. 
Keeping now a ^> fixed, and letting r vary from finite values to the deep UV 
regime, i.e. r = 0, the scale parameter r' governing the massless TBA systems 
will pass different regions. For the regime ln(2/r) < a/2 we see that the scaling 
functions (|36D and fl37f) are evaluated at r' > 1, whereas for ln(2/r) > a/2 they 
are taken at r' < 1. Thus, when performing the UV limit of the HSG model the 
massless TBA systems pass formally from the "infrared" to the "ultraviolet" regime 
with respect to the parameter r' . We emphasize that the scaling parameter r' has 
only a formal meaning and that the physical relevant limit we consider is still the 
UV limit r — > of the HSG model. However, proceeding this way has the advantage 
that we can treat the scaling function of the HSG model by the UV and IR central 
charges of the systems (B2T) and (K3T) as functions of r' 



[ CUV + CUV , 2 ^ 111 r 

Here we defined the quantities cjr '■= \vm r '-+oo c (r'), cjjv '■= lini r '-+o c Q (r') and 
cm, cuv analogously in terms of c (r'). 

In the case of cir + dm ^ cuv + cuv, we infer from ( p9|) that the scaling func- 
tion develops at least two plateaus at different heights. A similar phenomenon was 
previously observed for theories discussed in [T3], where infinitely many plateaus 



occurred which prompted to call them "staircase models" . As a difference, however, 
the TBA equations related to these models do not break parity. In the next sub- 
section we determine the central charges in (|39| ) by means of standard TBA central 
charge calculation, setting up the so-called constant TBA equations. 



4.2 Central charges from constant TBA equations 

In this subsection we will perform the limits r' — >• and r' — > oo for the massless 
systems ( |32| ) and (^) referring to them formally as "UV-" and "IR-limit" , respec- 
tively, keeping however in mind that both limits are still linked to the UV limit of 
the HSG model in the scale parameter r as discussed in the preceding subsection. 
We commence with the discussion of the "UV limit" r' — > for the subsystem 
(B^) . We then have three different rapidity regions in which the pseudo-energies are 
approximately given by 

( r'M a e ±e , for ±9 > -lnr' 

-£ 6 V\* for lnr'«0«-lnr' . (40) 

I -E b <P ab *L±(6), for ±0<lnr' 
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We have only kept here the dominant terms up to exponentially small corrections. 
We proceed analogously to the discussion as may be found in [0, |29|. We see that in 
the first region the particles become asymptotically free. For the other two regions 
the TBA equations can be solved by assuming the L-functions to be constant. 
Exploiting once more that the TBA kernels are centered at the origin and decay 
exponentially, we can similar as in section 3.1 take the L-functions outside of the 
integrals and end up with the sets of equations 



k-l 



JJ(1 + x±) Nab {l + xf) N '"b for lnr' < 9 < -lnr' (41) 

6=1 

k-l 

x a = Y[{l + x b f ab for ±0<lnr' (42) 

6=1 

for the constants x^ = e _e " ^ and x a = e _e " The N-matrices can be read off 
directly from the integral representations fl2"9| ) and (|30| ) 

N := — [<f> = l- 2(K su(k) )- 1 and N' := — [ if) = (K^)' 1 . (43) 
27T J 2n J 



Since the set of equations fl42|) has already been stated in the context of minimal 
ATFT and its solutions may be found in p5| , we only need to investigate the 
equations flU]) . These equations are simplified by the following observation. Sending 
9 to —9 the constant L-functions must obey the same constant TBA equation fl4"T| ) 
but with the role of and L~ interchanged. The difference in the masses mi, vn^ 
has no effect as long as ni\ ~ m2 since the on-shell energies are negligible in the 
central region lnr' <C 9 <C — lnr'. Thus, we deduce x+ = x~ =: x a and (|41| ) reduces 
to 



fc-i 

X n 



Y[(l + x b ) Nab with N = 1 - (fp«( fc ))-i . (44) 

6=1 

Remarkably, also these set of equations may be found in the literature in the context 
of the restricted solid-on-solid models P0fl . Specializing some of the general Weyl- 
character formulae in |2D| to the sw(3)fc-case a straightforward calculation leads 
to 

- (^) sin (Sgfl) Bi^(^l) 

x a = - ^ — ^ ^- - 1 and x a = - r - 1 . (45) 

^ m) sin sin (^) sin f^ff 



Having determined the solutions of the constant TBA equations fl4~T|) and ([P]) one 



can proceed via the standard TBA calculations along the lines of |2L [L3L E9[ and 
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compute the central charges from (|36|), (|37|) and express them in terms of Roger's 
dilogarithm function 



6 



k-l 



CUV 



lim c (r') = — 
6 



a=l 
fc-l 



2C 



X 



cuv = hm c (r' 



Using the non-trivial identities 



7T- 



a=l 



1 + X 

X a 



c 



X , , 



1 + Xa 



7T 



fc-1 



.1 /) 



1 + Xa 



+ X c 



k-l 



(46) 
(47) 



k + 3 



and 



TP 



.1 /; 



1 + x a 



k + 2 



(4E 



found in pi and |T3, we finally end up with 



cuv 



(fc-l)(4fc + 6) 
(fe + 3) (fc + 2) 



and 



cuv 



k + 2 



(49) 



For the reasons mentioned above cuv coincides with the effective central charge 
obtained from su(k) minimal ATFT describing parafermions |7j in the conformal 
limit. Notice that c uv corresponds to the coset (SU(3) k /U(l) 2 )/(SU(2) k /U(l)). 

The discussion of the infrared limit may be carried out completely analogous to 
the one performed for the UV limit. The only difference is that in case of the system 
(|3^) the constant TBA equations flHf ) drop out because in the central region the 
free energy terms becomes dominant when r' — > oo. Thus in the infrared regime 
the central charges of both systems coincide with &uv, 



c IR = lim Co(r') = c m = lim c G (r') = 2 



k + 2 



(50) 



In summary, collecting the results fl49|) and (50|), we can express equation (|39| ) 
explicitly in terms of the level k, 



c r, 




for 
for 



1< 2 < MS 
M|<2 



(51) 



We have replaced the limits in (^) by mass scales in order to exhibit the underlying 
physical picture. Here M~ is the smallest mass of an unstable bound state which 
may be formed in the process (a, i) + (b,j) — > (c,k) for Kf- ^ 0,2. We also used 

that the Breit-Wigner formula (|TTD implies that Mf ~ e "/ 2 for large positive a. 

First one should note that in the deep UV limit we obtain the same effective 
central charge as in section 3.1, albeit in a quite different manner. On the mathe- 
matical side this implies some non-trivial identities for Rogers dilogarithm and on 
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the physical ( |5T]) exhibits a more detailed behaviour than the analysis in section 3.1. 
In the first regime the lower limit indicates the onset of the lightest stable soliton 
in the two copies of complex sine-Gordon model. The unstable particles are on an 
energy scale much larger than the temperature of the system. Thus, the dynamical 
interaction between solitons of different colours is "frozen" and we end up with 
two copies of the SU(2)k/U(l) coset which do not interact with each other. As 
soon as the parameter r reaches the energy scale of the unstable solitons with mass 
Mg, the solitons of different colours start to interact, being now enabled to form 
bound states. This interaction breaks parity and forces the system to approach 
the SU(3)k/U(l) 2 coset model with central charge given by the formula in (H) for 
G = SU(3). 

The case when a tends to infinity is special and one needs to pay attention to 
the order in which the limits are taken, we have 

4 = lim lim c(r, a) 7^ lim lim c(r, a) = 6 . (52) 

k + 2 r-»0o--*ao cr^oor^O K + 3 

One might enforce an additional step in the scaling function by exploiting the 
fact that the mass ratio mi/iri2 is not fixed. So it may be chosen to be very large or 
very small. This amounts to decouple the TBA systems for solitons with different 
colour by shifting one system to the infrared with respect to the scale parameter 
r. The plateau then has an approximate width of ~ ln|m 1 /m 2 | (see figure 2). 
However, as soon as r becomes small enough the picture we discussed for m x ~ m 2 
is recovered. 



4.3 Restoring parity and eliminating the resonances 

In this subsection we are going to investigate the special limit a —>■ which is 
equivalent to choosing the vector couplings A± in (|T[) parallel or anti-parallel. For 
the classical theory it was pointed out in || that only then the equations of motion 
are parity invariant. Also the TBA-equations become parity invariant in the absence 
of the resonance shifts, albeit the S-matrix still violates it through the phase factors 
rj. Since in the UV regime a small difference in the masses mi and m 2 does not 
effect the outcome of the analysis, we can restrict ourselves to the special situation 
mi = m 2 , in which case we obtain two identical copies of the system 

*M + X>o6 + W * L b (e) = r M a cosh 6 . (53) 

6=1 

Then we have e o (0) = e { a\9) = e ( a\6), M a = M ( a ] = M { a ] and the scaling function 
is given by the expression 

c (r, a = 0) = ( ^J2 M * [ d0 L M cosh 6 ■ ( 54 ) 
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The factor two in comparison with ( |T9"D takes the two copies for % — 1,2 into account. 
The discussion of the high-energy limit follows the standard arguments similar to 
the one of the preceding subsection and as may be found in 0, p9"| . Instead of 
shifting by the resonance parameter a, one now shifts the TBA equations by lnr/2. 
The constant TBA equation which determines the UV central charge then just 
coincides with (f4ip. We therefore obtain 

12 fx \ k — 1 
lim lim c(r, a) — — \^ L ( — | = 6 . (55) 

7T 2 ^ \1+X a J k + 3 K ' 

Thus, again we recover the coset central charge (§) for G = SU(3), but this time 
without breaking parity in the TBA equations. This is in agreement with the results 
of section 3.1, which showed that we can obtain this limit for any finite value of a. 



4.4 Universal TBA equations and Y-systems 

Before we turn to the discussion of specific examples for definite values of the level 
k, we would like to comment that there exists an alternative formulation of the 
TBA equations (|T6|) in terms of a single integral kernel. This variant of the TBA 
equations is of particular advantage when one wants to discuss properties of the 
model and keep the level k generic. By means of the convolution theorem and the 
Fourier transforms of the TBA kernels <fi and ip, which can be read off directly from 
(p^) and (|30|), one derives the set of integral equations 

k-i 

e*(0) + Q k * U a {6 - = £ I ab Q k * ^ + V b ){e) . (56) 

6=1 

We recall that / denotes the incidence matrix of su(k) and the kernel f4 is found 
to be 

cosh(kt)/ 2) 



The on-shell energies have dropped out because of the crucial relation |pl | 

fc-i 



^/ afe M^ = 2cosf M l a , (5J 



6=1 

which is a property of the mass spectrum inherited from affine Toda field theory. 
Even though the explicit dependence on the scale parameter has been lost, it is 
recovered from the asymptotic condition 

ei(9)— * rMle ±e . (59) 
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The integral kernel present in (|56| ) has now a very simple form and the k dependence 
is easily read off. 

Closely related to the TBA equations in the form fl56|) are the following func- 
tional relations also referred to as Y-systems. Using complex continuation (see e.g. 



14|j for a similar computation) and defining the quantity Y^{6) = exp(— e l a {6)) the 



integral equations are replaced by 

k-i 

Yt(8 + %)Y*(6 ~ T ) = I 1 + Y H - **)] II t 1 + Yb(0rY Iab • ( 60 ) 

6=1 

The Y-functions are assumed to be well defined on the whole complex rapidity 
plane where they give rise to entire functions. These systems are useful in many 
aspects, for instance they may be exploited in order to establish periodicities in the 
Y-functions, which in turn can be used to provide approximate analytical solutions 
of the TBA-equations. The scaling function can be expanded in integer multiples 
of the period which is directly linked to the dimension of the perturbing operator. 

Noting that the asymptotic behaviour of the Y-functions is lim^oo Y*(6) ~ 
e -rM i a cosh 6^ we recover f or ^he Y-systems of the sw(fc)-minimal ATFT derived 



originally in In this case the Y-systems were shown to have a period related 
to the dimension of the perturbing operator (see fl8"3|)). We found some explicit 
periods for generic values of the resonance parameter a as we discuss in the next 
section for some concrete examples. 



5 Explicit examples 

In this section we support our analytical discussion with some numerical results and 
in particular justify various assumptions for which we have no rigorous analytical 
argument so far. We numerically iterate the TBA-equations ( |TED and have to choose 
specific values for the level k for this purpose. As we pointed out in the introduction, 
quantum integrability has only been established for the choice k > h. Since the 
perturbation is relevant also for smaller values of k and in addition the S-matrix 
makes perfect sense for these values of k, it will be interesting to see whether the 
TBA-analysis in the case of su(3)k will exhibit any qualitative differences for k < 3 
and k > 3. From our examples for the values k = 2, 3, 4 the answer to this question 
is that there is no apparent difference. For all cases we find the staircase pattern of 
the scaling function predicted in the preceding section as the values of a and x sweep 
through the different regimes. Besides presenting numerical plots we also discuss 
some peculiarities of the systems at hand. We provide the massless TBA equations 
(|32^) with their UV and IR central charges and state the Y-systems together with 
their periodicities. Finally, we also comment on the classical or weak coupling limit 
k — > oo. 
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5.1 The SU(3) 2 -HSG model 

This is the simplest model for the su(3)k case, since it contains only the two self- 
conjugate solitons (1,1) and (1,2). The formation of stable particles via fusing is 
not possible and the only non-trivial S-matrix elements are those between particles 
of different colour 

SH = SH = -1, S\l{9 -a) = -Sl\{9 + a)= tanh Ity- i\) . (61) 



Here we have chosen r] 12 = —rj 2 i = i- One easily convinces oneself that ( 51) satisfies 



indeed (|]) and (|1(]) . This scattering matrix may be related to various matrices which 
occurred before in the literature. First of all when performing the limit a — > oo 
the scattering involving different colours becomes free and the systems consists of 
two free fermions leading to the central charge c = 1. Taking instead the limit 
a — > the expressions in (|6l|) coincide precisely with a matrix which describes the 
scattering of massless "Goldstone fermions (Goldstinos)" discussed in |13|]. Apart 
from the factor i, the matrix Sfl(9)\ a=0 was also proposed to describe the scattering 
of a massive particle [|33| . Having only one colour available one is not able to set up 
the usual crossing and unitarity equations and in [Q the authors therefore resorted 
to the concept of "anti-crossing". As our analysis shows this may be consistently 
overcome by breaking the parity invariance. The TBA-analysis is summarized as 
follows 



unstable particle formation : c su(3) 2 = ~ = c uv + <ksv — ttt + 



6 .71 

_ = Cuv + Cuv — 1 — 

5 10 2 



no unstable particle formation : 2c su ( 2 ) 2 = 1 = cir + Cm = - + - . 

It is interesting to note that the flow from the tricritical Ising to the Ising model 



which was investigated in JT3[, emerges as a subsystem in the HSG-model in the 
form cuv ~^ c /i?- This suggests that we could alternatively also view the HSG- 
system as consisting out of a massive and a massless fermion, where the former is 
described by (PB]),(P2[) and the latter by ([J7|),([£|), respectively. 

Our numerical investigations of the model match the analytical discussion and 
justifies various assumptions in retrospect. Figure 1 exhibits various plots of the 
L-functions in the different regimes. We observe that for ln(2/r) < a/2, o ^ 
the solutions are symmetric in the rapidity variable, since the contribution of the 
■0 kernels responsible for parity violation is negligible. The solution displayed is 
just the free fermion L-function, L % {9) = ln(l + e~ rMl coshe '). Approaching more and 
more the ultraviolet regime, we observe that the solutions L l cease to be symmetric 
signaling the violation of parity invariance. The second plateau is then formed, 
which will extend beyond 6 = for the deep ultraviolet (see figure 1). The staircase 
pattern of the scaling function is displayed in figure 2 for the different cases discussed 
in the previous section. We observe always the value 6/5 in the deep ultraviolet 
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regime, but depending on the value of the resonance parameter and the mass ratio 
it may be reached sooner or later. The plateau at 1 corresponds to the situation 
when the unstable particles can not be formed yet and we only have two copies of 
sit (3) 2 which do not interact. Choosing the mass ratios in the two copies to be very 
different, we can also "switch them on" individually as the plateau at 1/2 indicates. 



1,0 - 

0,9 - 

0,8 - 

0,7 - 

0,6 - 

L (D (0) 0,5 - 

0,4 - 

0,3 - 
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Figure 1: Numerical solution for L^\9) of the su(3)2 related TBA-equations at different 
values of the scale parameter r and fixed resonance shift and mass ratio. 

The Y-systems (f)0|) for k = 2 read 

Yi(e + i^Yl(e-i^) = i + Y*(e-a ji ) 1.2. //./. (62) 



For cr = they coincide with the ones derived in |L3| for the "massless" subsystem. 
Shifting the arguments in (|62|) appropriately, the periodicity 



Yl{e + ^ + a^=Y({e) (63) 

is obtained after few manipulations. For a vanishing resonance parameter (|63|) 
coincides with the one obtained in g, [13|]. These periods may be exploited in a 
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series expansion of the scaling function in terms of the conformal dimension of the 
perturbing operator. 



1,2 
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Figure 2: Numerical plots of the scaling function for su(3)k, k = 2,3,4 as a function of 
the variable logr/2 at different values of the resonance shift and mass ratio. 

5.2 The SU(3) 3 -HSG model 

This model consists of two pairs of solitons (1, 1) = (2, 1) and (1, 2) = (2, 2). When 
the soliton (1, i) scatters with itself it may form (2, i) for % — 1, 2 as a bound state. 
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The two-particle S-matrix elements read 

Since soliton and anti-soliton of the same colour obey the same TBA equations we 
can exploit charge conjugation symmetry to identify e l {0) := e\(9) = e l 2 {6) leading 
to the reduced set of equations 

*{0) + i P *L i {0)-<p*ti{0-<T Si )=TM i <xxhO, ^) = - 1 4 ^ s S ||^ . (65) 
The corresponding scaling function therefore acquires a factor two, 

c(r, a) = ^ M% [ d0 cosh 9 U ( d ) ■ ( 66 ) 

This system exhibits remarkable symmetry properties. We consider first the situa- 
tion (7 = with mi = rri2 and note that the system becomes free in this case 

M (i) = M (2) = . M ^ £ (i)^ = £ (2)^ = rMcosh # . ( 67 ) 

meaning that the theory falls apart into four free fermions whose central charges 
add up to the expected coset central charge of 2. Also for unequal masses m\ ^ rri2 
the system develops towards the free fermion theory for high energies when the 
difference becomes negligible. This is also seen numerically. 

For o 7^ the two copies of the minimal A 2 -ATFT or equivalently the scaling 
Potts model start to interact. The outcome of the TBA-analysis in that case is 
summarized as 

6 4 

unstable particle formation : C W3) 3 = 2 = cuv + Qyy = — I — 

5 5 

8 4 4 

no unstable particle formation : 2c su ( 2 ) 3 = — = cm + cm — — I — . 

5 5 5 

As discussed in the previous case for k = 2 the L-functions develop an additional 
plateau after passing the point ln(2/r) = a/2. This plateau lies at In 2 which is the 
free fermion value signaling that the system contains a free fermion contribution 
in the UV limit as soon as the interaction between the solitons of different colours 
becomes relevant. Figure 2 exhibits the same behaviour as the previous case, we 
clearly observe the plateau at 8/5 corresponding to the two non-interacting copies 
of the minimal A 2 -ATFT. As soon as the energy scale of the unstable particles is 
reached the scaling function approaches the correct value of 2. 
The Y-systems ( pO]) for k = 3 read 

■ / 7T\ • / 7T\ 1 + Y, j J9 + an) 

{o + ij) y;, (» - « s ) = w iJ = h 2 ' '* j ' (68) 
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Once again we may derive a periodicity 

Y^{e + 2m + a ji )=Y^{6) (69) 
by making the suitable shifts in ( |5"8"D and subsequent iteration. 

5.3 The SU(3) 4 -HSG model 

This model involves 6 solitons, two of which are self-conjugate (2,1) = (2,1), 
(2,2) = (2,2) and two conjugate pairs (1,1) = (3,1), (1,2) = (3,2). The cor- 
responding two-particle S-matrix elements are obtained from the general formulae 
(§ and (|) 

/ (2), (3) e (l) e -(2) e \ 
S u (6) = (3)e(l) 6 (2)1 (3),(1), (70) 

V -(2), (3),(1), (2), J 

for soliton-soliton scattering with the same colour values and 

/ Vij(-l)e ^{-2)6 i&(-3), \ 
SF(e -a ij )=\ r]l(-2)e -(-3)g(-l) e i&(-2), (71) 

V 4(- 3 )e vU~2)e / 

for the scattering of solitons of different colours with ry 12 = e*5. In this case the 
numerics becomes more involved but for the special case mi = m<i one can reduce 
the set of six coupled integral equations to only two by exploiting the symmetry 
l£\6) = La (—9) and using charge conjugation symmetry, L\{9) = L l 3 (8). The 
numerical outcomes, shown in figure 2 again match, with the analytic expectations 
(p|) and yield for ln(2/r) > cr/2 the coset central charge of 18/7. In summary we 
obtain 

unstable particle formation : c su ( 3 ) 4 = — = C\jy + Cuv — -y + 1 

no unstable particle formation : 2c SM ( 2 ) 4 = 2 = Cj# + Cm = 1 + 1 , 

which matches precisely the numerical outcome in figure 2, with the same physical 
interpretation as already provided in the previous two subsections. 

5.4 The semi-classical limit k — > oo 

As last example we carry out the limit k — > oo, which is of special physical interest 
since it may be identified with the weak coupling or equivalent ly the classical limit, 
as is seen from the relation h(3 2 = l/k + 0(l/k 2 ). To illustrate this equivalence 
we have temporarily re-introduced Planck's constant. It is clear from the TBA- 
equations that this limit may not be taken in a straightforward manner. However, 
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we can take it in two steps, first for the on-shell energies and the kernels and finally 
for the sum over all particle contributions. The on-shell energies are easily computed 
by noting that the mass spectrum becomes equally spaced for k — ► oo 



Ml = M, 



k—a 



rrii it a 
—2 sm — a rrii 
ir(3 2 k 



k 

a < — . 
2 



(72) 



For the TBA-kernels the limit may also be taken easily from their integral repre- 
sentations 



U*) -— 2tt*(0) [5 ab -2[K ab 



su(k) 



-1 



and ^{0)^2*6(0) [K s 2 {k) 



k— >oo 



(73) 

when employing the usual integral representation of the delta-function. Inserting 
these quantities into the TBA-equations yields 



fc-i 



4(0) 



r am. 



6=1 



^ab 



fc-1 



Hi 



6=1 



(74) 

We now have to solve these equations for the pseudo-energies. In principle we could 
proceed in the same way as in the case for finite k by doing the appropriate shifts 
in the rapidity. However, we will be content here to discuss the cases a —>■ and 
a — > oo, which as follows from our previous discussion correspond to the situation 
of restored parity invariance and two non-interacting copies of the minimal ATFT, 
respectively. The related constant TBA-equations fl4"2|) and (041) become 



oo : x a 



(a + 1) 2 
fc^oo a(a + 2) 



1 and a — > : x n 



(o + l)(q + 2) 
k-^oo a(a + 3) 



1. (75) 



The other information we may exploit about the solutions of (|f4|) is that for large 
rapidities they tend asymptotically to the free solution, meaning that 



o — > 0, oo : 



+±oo 



ln(l + e" 



-r arm cosh 8\ 



(76) 



We are left with the task to seek functions which interpolate between the properties 
(f75p and (|76|) . Inspired by the analysis in we take these functions to be 



a 



a 



■yz 



: 



L l M = In 
V a {9) = In 



sinh 2 (^ii raij cosh 9) 



sinh (|- nrij cosh 9) sinh rrrij cosh 9\ 

sinh (^p rnij cosh 9) sinh (^i 2 - rm^ cosh #) 
sinh (| rmj cosh 9) sinh f^^r rm; cosh 9) 



(77) 
(78) 



24 



The expression ([77]) coincides with the expressions discussed in the context of the 



breather spectrum of the sine-Gordon model and (|78|) is constructed in analogy. 



We are now equipped to compute the scaling function in the limit k — > oo 

o 2 fc-1 

c(r, a) = lim 4 E / ^ cosh 61 E M * L ^) " ( 79 ) 

i=l ^ a=l 

Using (^), (f77|) and (|78|) the sum over the main quantum number may be computed 



directly by expanding the logarithm. We obtain for k — > oo 
C (r)| CT=oc = ^^/"^m J cosh^ln(l-e- rm!COsh9 ) (80) 

1 i=:r 



-4^^/'^m i cosh0[ln(l-e- rm!COsh9 ) + ln(l - e - r2m * coshe )]. (81) 
^ i=:r 



c(r) CT=0 



Here we have acquired an additional factor of 2, resulting from the identification of 



particles and anti-particles which is needed when one linearizes the masses in (|72l) . 
Taking now the limit r —>■ we obtain 

no unstable particle formation : 2 c^^)^ = 4 (82) 
unstable particle formation : c sn ( 3 ) tx) = 6 . (83) 



The results (|30"D, (|82"D and (|3lD, Q83|) allow a nice physical interpretation. We 
notice that for the case cr — > oo we obtain four times the scaling function of a 
free boson. This means in the classical limit we obtain twice the contribution of 
the non-interacting copies of SU(2) 00 /U(1), whose particle content reduces to two 
free bosons each of them contributing 1 to the effective central charge which is in 
agreement with (H). For the case a —> we obtain the same contribution, but in 
addition the one from the unstable particles, which are two free bosons of mass 2mj. 
This is also in agreement with (||). 

Finally it is interesting to observe that when taking the resonance poles to be 
9r = a — in jk the semi-classical limit taken in the Breit-Wigner formula (|TT|) leads 
to m| = (raj + rrij) 2 . On the other hand (pi|) seems to suggest that m~ k = 2rrii, 
which implies that the mass scales should be the same. However, since our analysis 
is mainly based on exploiting the asymptotics we have to be cautious about this 
conclusion. 



6 Conclusions 

Our main conclusion is that the TBA-analysis indeed confirms the consistency of 
the scattering matrix proposed in 0. In the deep ultraviolet limit we recover the 
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Gk/U(l) -coset central charge for any value of the 2£ — 1 free parameters entering 
the S-matrix, including the choice when the resonance parameters vanish and parity 
invariance is restored on the level of the TBA-equations. This is in contrast to the 
properties of the S-matrix, which is still not parity invariant due to the occurrence 
of the phase factors 77, which are required to close the bootstrap equations Q. 
However, they do not contribute to our TBA-analysis, which means that so far 
we can not make any definite statement concerning the necessity of the parity 
breaking, since the same value for the central charge is recovered irrespective of the 
value of the ex's. The underlying physical behaviour is, however, quite different as 
our numerical analysis demonstrates. For vanishing resonance parameter the deep 
ultraviolet coset central charge is reached straight away, whereas for non-trivial 
resonance parameter one passes the different regions in the energy scale. Also the 
choice of different mass scales leads to a theory with a different physical content, 
but still possessing the same central charge. To settle this issue, it would therefore 
be highly desirable to carry out the series expansion of the scaling function in r and 
determine the dimension A of the perturbing operator. It will be useful for this to 
know the periodicities of the Y-functions. We conjecture that they will be 

Y: (9 + ztt(1 - A)- 1 + an) = Yi{6), (84) 

which is confirmed by our sii(iV)-examples. For vanishing resonance parameter and 
the choice g = su(2), this behaviour coincides with the one obtained in f28fl . This 
means the form in (Q) is of a very universal nature beyond the models discussed 
here. 

We also observe from our sii (A)-example that the different regions, i.e. k > h 
and k < h, for which quantum integrability was shown and for which not, respec- 
tively, do not show up in our analysis. 

It would be very interesting to extend the case-by-case analysis of section 5 to 
other algebras. The first challenge in these cases is to incorporate the different 
resonance parameters. 
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